Deconfinement Phase Transition in a 3D Nonlocal U(l) Lattice Gauge Theory 
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We introduce a 3D compact U(l) lattice gauge theory having nonlocal interactions in the temporal 
direction, and study its phase structure. The model is relevant for the compact QED3 and strongly 
correlated electron systems like the t-J model of cuprates. For a power-law decaying long-range 
interaction, which simulates the effect of gapless matter fields, a second-order phase transition takes 
place separating the confinement and deconfinement phases. For an exponentially decaying inter- 
action simulating matter fields with gaps, the system exhibits no signals of a second-order transition. 
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The U(l) lattice gauge theory (LGT) in three dimen- 
sions (3D) coupled to matter fields describes various in- 
teresting physical systems. The compact QED3 is just a 
such system and its phase structure has been studied by 
various methods 1 . In condensed matter physics, interest- 
ing observations were made that the strongly-correlated 
electron systems in two dimensions like the antiferromag- 
netic Heisenberg spin model, the t-J model of high-T c 
cuprates, and the fractional quantum Hall states are de- 
scribed by 3D U(l) gauge theories due to the introduction 
of auxiliary collective fields 2-5 . 

For the case without matter-field couplings, Polyakov 6 
showed that the 3D compact U(l) gauge theory is always 
in the confinement phase due to the monopole (instan- 
ton) condensation. For the case with couplings to mat- 
ter fields, there is still no consensus on the question of 
whether the system in two spatial dimensions exhibits a 
phase transition into a deconfinement phase. Probably 
the answer depends on the properties of coupled mat- 
ter fields. This question is important because the "frac- 
tionalization" of electrons may be interpreted as the de- 
confinement phenomenon of U(l) gauge dynamics. For 
the t-J model, the possibility of charge-spin separation 
(CSS) is of great interest since it may explain various 
anormalous behaviors of the metallic state of cuprates . 
In Ref. 8 it was argued that the CSS takes place below 
certain critical temperature (T) as a deconfinement (per- 
turbative) phase of an effective U(l) LGT which is de- 
rived by the hopping expansion of spinons and holons in 
the slave-particle representation at finite T. In related 
works, Nagaosa 9 argued that the deconfinement phase is 
possible above some T, whereas Nayak 10 argued that de- 
confinement does not occur at any T. The deconfinement 
phase at T = for systems with gapless excitations are 
supported in Ref. 11 and denied in Rcf. 12 . 

In this Letter, we introduce and study a LGT with non- 
local interactions in order to investigate the phase struc- 
ture of compact U(l) gauge theories coupled to matter 



PACS numbers: 11. 15. Ha, 12.38.Gc, 71.27. +a 

fields on the 3D lattice (2D system at T = 0). We first 
consider the cases of massless and massive relativisitic 
matter fields. Then we apply the model to the nonrel- 
ativistic electron systems. The results of this paper for 
gapless excitations in electron systems shall complement 
our previous result of CSS 8 because the hopping expan- 
sion employed there may be inadequate for massless ex- 
citations at T = 0. 

Let us start with the path-integral representation of 
the partition funciton Z of gauge field U xfJl and (bosonic 
and/or fermionic) matter fields <p x , 



X X . fJ, 

A = - ^ 4> x T xy {U)(j) v + A v , 



Au = q^2 (UxMx+i>^U x+ f l>u U xtl + c.c.), (1) 

where x = (xq, X\, X2) is the site-index of the 3D lat- 
tice of the size V — NqN%N2 with the periodic bound- 
ary condition, fi(= 0,1,2) is the (imaginary) time and 
spatial direction index, <p x is the matter field on x, 
U XI1 = exp(i0 aM ) (-7T < 6 XIM < 7r) is the U(l) gauge vari- 
able on the link (x,x + /t), T xy (U) represents the local 
couplings of (f> x to U xfi . 

By integrating over <p x , we obtain an effective gauge 
theory, 



Z 



J\dU x 



. exp 



/Tr log T xy (U)+Au 



(2) 



where / is a parameter counting the statistics and inter- 
nal degrees of freedom of <p x . Due to the (Tr log) term, 
the effective gauge theory becomes nonlocal. For rela- 
tivistic matter fields, it is expanded as a sum over all the 
closed random walks 1Z (loops including backtrackings) 



1 



on the 3D lattice which represent world lines of particles 
and antiparticles as 



Tr log T xy (U) = J2 



7 



L[R] 



L[K] 



n 



(3) 



(i/i)GR 



L[R] is the length of TZ, and 7 = (6 + m 2 ) -1 is the hop- 
ping parameter (6 is the number of links emanating from 
each site, and m is the mass of the matter field in unit of 
the lattice spacing). For the constant gauge-field config- 
uration Uxn = 1, the expansion in (3) is logarithmically 
divergent ~ logm as m — > due to the lowest-energy 
zero-momentum mode. 

Below we shall study a slightly more tractable model 
than that given by Eq.(2). It is suggested from the formal 
hopping expansion (3), and obtained by retaining only 
the rectangular loops extending in the temporal direction 
in the loop sum and choosing their coefficients optimally 
as follows; 



Z T = f Y[dU x ^exp(A T ), 



2 N 

A T = CrO^vr + V *,i,r) + A S, 

x 2—1 r — 1 
r-1 

V Xj i tT = U x+r Q i JJ [U x+k QfiU x+ ~ i+k Q }U x i, 
k=0 

As = g\^{Ux2U x+ % 1 U x+l2 U xl + c.c). 



(4) 



V x> i >r is the product of U xfl along the rectangular (x, x + 
i, x + i + rO, x + rO) of size (1 x t) in the (i — 0) plane. 
In As, we have retained only the nearest-neighbor spa- 
tial coupling. For the nonlocal coupling constant c T , we 
consider the following two cases; 



power — law decay (PD), 
exponential decay (ED). 



(5) 



The power —1 in the PD case in (5) reflects the ef- 
fect of massless excitations. In fact, this c T gener- 
ates a logarithmically divergent action for £/ XjU — 1 ex- 
plained below Eq.(3) as one can see from the relation, 
^ exp(— mr)T _1 ~ log(l/m). The action for m = is 
then proportional to X)r T 1 — log-^o for finite No. On 
the other hand, the ED model contains m and simulates 
the case of massive matter fields 13 . 

We made Monte Carlo simulations to determine the 
phase structure of this model. We consider the isotoropic 
lattice, N^ = N, with the periodic boundary condition 
up to N = 32, where the limit N — > 00 corresponds to 
the system on a 2D spatial lattice at T — 0. For the mass 
of the ED model, we set m = 1. For the spatial coupling 
A scaled by g, we consider the two typical cases A = 
(i.e., no spatial coupling) and A = 1. 

First, we calculate the following average E ("internal 
energy") and the fluctuation C ("specific heat") of Aj-; 
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FIG. 1. Internal energy E 
the action with A = 1 vs 
N = 8(f), 16(*), 24(B); (a,c) PD model, (b,d) ED model. 
The solid lines in (a) and (b) are the large-;? expansion. In the 
PD model, strong iV dependence is observed in E at large g 
and in the developing peak of C. They indicate a second-order 
phase transition in the PD model. 



E = -{A T )/V, C=((A T -(A T )) 2 )/V. 



(6) 



For small g, the high-temperature expansion (HTE) gives 
Z ~ exp[ 5 2 (2Q 2 + A 2 )V r ] (Q k = £ T <*), whereas for large 
g, the low-temperature expansion (LTE) around V X ^ :T = 
1 gives Z ~ exp[(4 ff Q 1 + 2.gA - \ogg)V]. 

In Fig.l, we present E,C for A = 1 vs the nonlocal 
coupling g. They agree with the above HTE and LTE. 
In the PD model, E of Fig. 1(a) connects the HTE result 
and LTE result, showing that V X i, r ~ 1 for large g. C of 
Fig. 1(c) shows that its peak develops as N increases. The 
finite-size scaling analysis shows C of Fig. 1(c) fits well in 
the form, C(g,N) = N°/ u <t){N x / v e),e= (g-g c )/g c with 
v = 1.2 ~ 1.3, a/v = 0.25 - 0.26, g c = 0.10 - 0.12. 
These results indicate that the PD model exhibits a 
second-order phase transition separating the disordered 
(confinement) phase and the ordered (deconfinement) 
phase at g = g c . This transition will be confirmed later 
by the measurement of Polyakov lines. On the other 
hand, the peak of C in the ED model does not develop as 
N increases, showing no signals of a second-order tran- 
sition. It may have a higher-order transition or just a 
crossover. Simulations of the models with A = give 
similar behaviors of E, C, preserving the above phase 
structure for A = 1. 

To study the nature of gauge dynamics, we calculated 
the spatial correlations of Polyakov lines P x± ; 



JV ( , 



Px ± = Y[ U x±,x ,0, (x± = {x 1 ,X 2 )), 

f P (x ± ) = (P x± P Q ). 



(7) 



Since the present model (4) contains no long-range inter- 
actions in the spatial directions, fp(x±) is expected 
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FIG. 2. Correlations of Polyakov lines, fp(x±), vs 
(a) PD (A = 0, N=16) with g = 0.4,0.3,0.2,0.1,0.02 from 
above, (b) ED (A = 0, N=16) with g = 2.5,2.0, 1.5, 1.0,0.5 
from above. (c) and (d) show the order parameter 
P = [/p(z± AX )] 1/2 vs g for the PD model; (c) A = and 
(d) A = 1. They exhibit a long-range order for g > g c ~ 0.15 
in the PD model. 



to supply us with a good order parameter to detect a 
possible confinement-deconfinement transition. The de- 
confinement phase is characterized by small fluctuations 
of U xfJl and therefore by an order in fp(xx). 

In Fig. 2, we present fp(x±) with A = 0. The PD 
model of Fig. 2(a) clearly exhibits an off-diagonal long- 
range order, i.e., lim^^oo fp(x±) ^ for g > 0.20, 
whereas the ED model of Fig. 2(b) does not for all g's. 
To see this explicitly, we plot in Fig.2(c,d) the order pa- 
rameter p = [/p(x r ^ AX )] 1 / 2 for the PD model, where 
x max = jyj ^2 i s the distance at which fp becomes min- 
imum due to the periodic boundary condition, p starts 
to develop continuously from zero at g — g c ~ 0.15. The 
size dependence of p shows a typical behavior of a second- 
order transition. Thus the gauge dynamics of the PD 
model is realized in the deconfinement phase for g > g c , 
whereas it is in the confinement phase for g < g c . In 
contrast, the ED model stays always in the confinement 
phase. These results including the value of g c are consis- 
tent with those derived from the data of E, C in Fig.l. 

To see the details of gauge dynamics, we measured 
the instanton density p{x), an index for disorder of 
U xfJ ,. We employ the definition of p(x) in U(l) LGT 
by DeGrand and Toussaint 14 . For the local 3D com- 
pact U(l) LGT without matter fields, the average density 
p = (J2 X \p( x )\)/V was measured in Ref. 15 . In Fig. 3 we 
present p vs g. It decreases as g increases more rapidly 
in the PD model than in the ED model. This behavior 
of p is consistent with the result that the PD model ex- 
hibits a second-order transition, while the ED model does 
not. The A coupling enhances the rate of decrease in p as 
one expects since the spatial coupling favors the ordered 
deconfinement phase. In the ED model with A = 1, p 
is fitted by oc exp(— eg) with a constant c in the dilute 
(large g) region, and the smooth increase for smaller g 
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FIG. 3. Average instanton density p vs g. (a)PD(A 
N = 8), (b)PD(A = 0, N = 16), (c)PD(A = 1, N = 8), 
(d)PD(A = 1, N = 16), (e)ED(A = 0, iV = 8,16), 
(f)ED(A = 1, N = 8,16). The solid curve oc exp(-cg) fits 
(f) at large g. 




FIG. 4. Snapshots of instanton configuration p(x) for 
N=16. (a) PD(A = l,g = 0.15) and (b) PD(A = l,g = 0.30). 
The light cubes for p(x) — 1 and the dark cubes for p(x) = — 1. 



indicates a crossover from the dilute gas of instantons to 
the dense gas, just the behavior similar to the case of 
pure and local LGT 6 ' 15 . 

In Fig. 4 we show snapshots of p(x) for the PD model 
with A = 1. Fig. 4(a) is a dense gas and Fig. 4(b) is a dilute 
gas. They are separated at g c ~ 0.20, the location of the 
peak of C for N = 16. In Fig. 4(b), instantons mostly 
appear in dipole pairs at nearest-neighbor sites, p(x) = 
1, p(x ± p) = — 1, while in Fig. 4(a), they appear densely 
and it is hard to determine their partners. In both cases, 
the distributions p(x) have no apparent anisotoropies like 
column structures. However, the orientations of dipoles 
in Fig. 4(b) are mostly (~ 92%) in the temporal direction 
as expected from Eq.(4). 

For ordinary pure and local gauge systems, Wilson loop 
W[C] = (Y\ c U xt _,) along a closed loop C on the lattice is 
used as an order parameter to study the gauge dynamics; 
W[C] obeys the area law in the confinement phase and 
the perimeter law in the deconfinement phase; 



W[C] 



exp(—aS[C]), area law, 
exp(— a'L[C]), perimeter law, 



(8) 



where S[C] is the minimum area of a surface, the bound- 
ary of which is C. For a local LGT containing matter 
fields, W[C] cannot be an order parameter because the 



3 



matter fields generate the terms Yic with coeffcients 
~ exp(— 6L[C]) in the effective action. However, in the 
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FIG. 5. Wilson loops (N — 32) in the 1-2 plane at large g vs 

L[C] or S[C]. (a) PD(A = l,g = 0.25), (b) ED(A = l,g = 1.5). 

The PD model seems to prefer the perimeter law, whereas the 

ED model prefers the area law. 



present model (4) , the nonlocal terms are restricted only 
along the temporal direction, so it is interesting to mea- 
sure W[C] for the loops lying in the spatial (1-2) plane. 

In Fig. 5, we plot W[Cj. For the PD model in Fig. 5(a), 
the data at g = 0.25 seem to prefer the perimeter law. 
For the ED model in Fig. 5(b), the area law fits W[C] 
better than the perimeter law at g = 1.5; a considerably 
larger value than g ~ 1.0 at the peak of C. This suggests 
the area law holds in the ED model at all g. These ob- 
servations are consistent with the previous results on the 
(non)existence of a phase transition. Wilson loops in the 
spatial plane are useful to study the gauge dynamics of 
the present model. 

We have observed that the nonlocal couplings along the 
temporal direction in the PD model have sufficient effect 
of suppressing fluctuations of U xfl to produce the decon- 
finement phase. This result strongly suggests a decon- 
finement transition in the original model (2) with mass- 
less matter fields, because the isotropically distributed 
nonlocal couplings of the original model should have sim- 
ilar effect. In such a possible deconfinement phase of the 
original model, perturbation theory may be applicable, 
which predicts the potential energy between two charges 
as V(r) ~ r _1 , a weaker one than the 3D Coulomb po- 
tential V{r) ~ log(r). 

Let us turn to the nonrelativistic case. For the t-J 
model, by using the hopping expansion of holons and 
spinons at finite T with the continuous imaginary time, 
we derived an effective gauge theory, which is highly non- 
local in the temporal direction. The obtained effective 
theory has a similar action as Eq.(4) with c T =constant 
and g n where n is the density of matter fields (holons 
and spinons) 8 . The nonlocal correlations of the effective 
gauge field like c T ^constant come from the fact that the 
nonrelativistic fermions contain a higher density of low- 
lying excitations compared to Dirac fermions, i.e., the 
existence of the Fermi surface (or line). Although the 
above effective gauge model is obtained for the system at 
finite T, we expect that a similar gauge model appears 



as an effective model at T = 0. Then it is interesting and 
also important to investigate the gauge model (4) with 
c r =constant, the nondecaying (ND) model. The ND 
model should shed some light on the anormalous normal 
state of cuprates because the Fermi line exists there. 

We also made a Monte Caro simulation of the ND 
model, and obtained a phase transition similar to the 
PD model. However, the developing peak of C shifts 
to smaller g as N increases more quickly than the PD 
model as g c ~ 0.1 (JV = 8),0.045(7V = 16),0.03(JV = 24) 
for A = 1. It seems likely that g c — ► as N — ► oo, that is, 
the deconfinement phase dominates for all g(> 0). This 
may be related with the diverging coefficient Q2(oc N) 
in HTE, which implies that the radius of convergence is 
zero. This dominance of deconfienment phase of the ND 
model at T = may support the CSS at finite T, which is 
consistent with the result of Ref. 8 . In contrast to the ND 
model, the PD model has a finite limit of Q%{= 1.645), 
and has a finite region < g < g c of the confinement 
phase. 
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